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Целью данной работы является исследование актуальной за- 
дачи коммивояжера, которая является МР сложной задачей 
дискретной оптимизации. Показано, что для достижения по- 
ставленной цели пригодны лишь эвристические методы. Для 
решения задачи представлен результат совместного использо- 
вания муравьиного (МА) и генетического (ГА) алгоритмов. 
Особенностью предложенного генетического алгоритма явля- 
ется то, что задача решается только с помощью различных 
мутаций (без кроссовера). Исследуемый ГА был улучшен вве- 
дением элитарной стратегии. Испытания проводились на гра- 
фах средней и большой размерности. Элитная особь, получае- 
мая с помощью муравьиного алгоритма, в среднем улучша- 
лась на 11%. Показано, что эффективность генетического 
алгоритма напрямую зависит от количества особей в поколе- 
нии и количества итераций алгоритма. Ввод элитарной страте- 
гии улучшил получаемые значения целевой функции более 
чем в 2 раза. Увеличение числа запусков МА при выборе эли- 
ты позволило повысить эффективность алгоритма приблизи- 
тельно на 2 %. 
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Введение. Задача коммивояжера (ЗК) является МР сложной задачей дискретной оптимизации и определяется 


следующим образом: для заданного полного взвешенного графа С = (У, Е, О) с множеством вершин У = {у,}„, множе- 


ством ребер Е и матрицей весов О = {4,;;} х, необходимо найти Гамильтонов цикл (у, Ур,...,Уш), где (1,... 12) — пере- 


становка на множестве {1,2.,...,п 1}, ау; с У, с минимальной суммой весов дуг [1]: 


п-1 


уа.. 


Л! 


111 


+4, —> шт 


Вершины графа С часто называются городами, а любой Гамильтонов цикл называется маршрутом коммивоя- 
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жера. Очевидным решением задачи является метод полного перебора. Выйдя из первого города, коммивояжер может 
направиться в один из (и-1) городов, откуда в (п-2) оставшихся городов и т. д., пока не останется один- 
единственный город, откуда он вышел. Таким образом, общее число маршрутов, подлежащих просмотру, составляет 
(п-1):(и-2)...2:1= (п- 1) всевозможных вариантов. С увеличением количества городов это значение быстро воз- 
растает и уже при 15 городах достигает астрономических цифр. Поиск точных и приближенных способов решения 
задачи о коммивояжере остается актуальным и с теоретической, и с практической точки зрения, так как не найдено, а 
возможно, и не существует точных алгоритмов, решающих ЗК за короткое (полиномиальное) время. Поэтому для ре- 
шения задачи коммивояжера эффективнее использовать не точные, а эвристические методы. Применение эвристиче- 
ских алгоритмов к какой-либо практической задаче обычно позволяет сформулировать рекомендации, которые намно- 


го лучше произвольного решения и которые наиболее близки к лучшему варианту. 


Обзор алгоритмов решения. В данной статье решение задачи коммивояжера произведено модифицирован- 
ным генетическим алгоритмом (ГА) с элитной особью. Начальное значение элиты получается при использовании п 
итераций муравьиного алгоритма (МА), который также является приближенным и существует в рамках научного 
направления Мабига! Сотрийп® — «Природные вычисления», которое объединяет математические методы и механиз- 
мы, обеспечивающие многолетнюю адаптацию биоценозов к окружающей среде [2—5]. Использованный МА также 
улучшен применением элитной стратегии. Кроме того, в данной версии МА используется модификация — «правило 
псевдослучайных пропорций», которая поддерживает баланс между использованием накопленных знаний и исследо- 
ванием новых решений. 

Рассматриваемый муравьиный алгоритм предполагает выполнение следующих шагов [6]. 

1. Инициализация: создание муравьев (каждый муравей помещается в начальную точку). Задается начальный 
уровень феромона (обычно небольшое положительное число) для того, чтобы на начальном этапе вероятности перехо- 
да в вершины не принимали нулевые значения. 

2. На втором шаге определяются вероятности перехода К-го муравья из вершины 1 в вершину ]. Для этого ис- 
пользуем следующую формулу (1): 

Е а 
ь] 
ро (1) 
У -т. — 
о ь/ 
р, 
где г; (0) — уровень феромона; [,, — расстояние между вершинами 1 и ]; Виа — константные параметры. 

При а = 0 выбор ближайшего города наиболее вероятен, то есть алгоритм становится жадным. При В = 0 вы- 

бор происходит только на основании феромона. 


3. На третьем шаге алгоритма применяется правило псевдослучайных пропорций, которое определяется про- 


стой формулой (2) и зависит от параметра 4: 


пах ( [ти }, 954 
р 7еаЙомеа, (2) 


Г, 42 40 


где 4 — случайное число, распределенное по нормальному закону на отрезке [0,1]; 40 — параметр баланса между ис- 
пользованием накопленных знаний и исследованием новых решений (0 < 40 < 1); г — случайный город, выбранный на 
основе вероятностей, расчет которых был произведен по формуле (1). 

4. Пройдя ребро (1, 7), муравей откладывает на нем некоторое количество феромона, которое должно быть свя- 


зано с оптимальностью сделанного выбора: 
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Здесь О — параметр, имеющий значение порядка цены оптимального решения; 7, (1) — маршрут, пройден- 


ный муравьем К к моменту времени (; [^(Р) — длина этого маршрута. 
5. Дополнительным улучшением алгоритма является ввод элитных муравьев. После каждой {-й итерации ал- 
горитма происходит дополнительное обогащение феромоном ребер, входящих в маршрут лучшего на данный момент 


муравья. 


е/ [2% (#),если(е, 7) ЕТ” 


0, если(1, 1) ® 1. 


Дт, — 


Особенностью генетического алгоритма, приведенного в данной статье, является то, что оптимальный марш- 
рут находится посредством модели Гольдберга с использованием только лишь модифицированных операторов мута- 
ции (ОМ) без применения кроссовера. Мутация — случайное изменение одного или нескольких генов в хромосоме [7]. 
Оператор мутации позволяет разнообразить генотип в поколении, тем самым исключая возможность «застревания» в 
ловушке локальных минимумов и повышая точность работы алгоритма. Вероятность мутации р„ является фиксиро- 
ванным случайным числом на отрезке [0; 1]. Вероятность мутации в предложенном алгоритме составляет 100 %. В 
ходе вычислительных экспериментов, изложенных в работе [8], было определено, что наиболее близкие к оптимуму 
длины маршрутов получаются при использовании простого несимметричного двухточечного ОМ. 

Несимметричный двухточечный ОМ выполняется следующим образом. 

1. В хромосоме Р = (а, а», а, ... ‚аи, авт, а.) определяются случайным образом две позиции (например, а› и а, )). 

2. Гены, соответствующие выбранным позициям, переставляются, и формируется новая хромосома 
Р'= (а, ат, а, ... ‚, ак 2, а2, а»). 

Для данной работы была выбрана стратегия «Сравнение мутирующей особи с предком, а далее — лучшей из 
них со случайной особью в поколении», которая показала себя как одна из лучших стратегий [9-10]. 

Алгоритм реализуется следующим образом. 

1. Получаем значение элитной особи при помощи п запусков МА. 

2. Ставим элитную особь на первое место в поколении. 

3. Если т — номер особи, подверженной мутации, то т = 2. 

4. Применяем к т-й особи оператор мутации. 

5. Вычисляем значение целевой функции для мутирующей т-й особи и сравниваем его со значением целевой 
функции предка, а затем лучшее из них сравниваем со значением целевой функции случайно выбранной особи из по- 
коления. Особь с наилучшим значением целевой функции добавляется в следующее поколение. 

6. Устанавливаем т = т + | и повторяем шаг 5, пока каждая особь из поколения не подвергнется мутации. 

7. Определяем особь с наилучшей приспособленностью (значением целевой функции) в поколении [и и ста- 
вим ее на первое место. 

8. Проверяем, сколько итераций подряд значение /„„ было неизменно. Если /„„ было неизменным Ма ите- 
раций подряд, то конец работы, если меньше, чем №, то переходим к шагу 3. 


Аналитически сделать ВЫВОДЫ об эффективности данного алгоритма нельзя, поэтому был поставлен вычисли- 


тельный эксперимент. Параметры использованного МА: число особей на итерации — 50, одна элитная особь, @=1, 


В =1. Результаты вычислительного эксперименты для графа пВае229 (29 вершин, известное значение оптимума 1610) 


представлены в табл. 1, а для графа — Бга7158 (58 вершин, известное значение оптимума 25395) при выборе элитной 


особи из трех запусков МА — в табл. 2, и при выборе из пяти запусков МА — в табл. 3. 
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Таблица 1 


Результаты вычислительного эксперимента граф — пВаез29 (число опытов — 35) 


Число Выбор элитной особи из 3 Выбор элитной особи из 5 


повторов/ Луч Сред. Сред. знач. Сред. % Луч Сред. Сред. знач. Сред. % 
особей шее знач. элиты улучш. шее знач. элиты улучш. 


1627 1782 1949 10.41 1615 1756 1921 10.25 
1620 | 1746 1926 11,18 1610 | 1751 1901 
1610 | 1709 1867 1610 | 1757 1924 10,32 


Из табл. 1 видно, что элитная особь в среднем улучшается на = 10 %. Лучшее из полученных значений отстает 





от известного оптимального в худшем случае всего на = 1 %. А в большинстве случаев достигает оптимума. Среднее 
значение отстает от оптимума всего на 8,6 %. 


Таблица 2 
Выбор элитной особи из трех (граф — га71158) 


Число 20 запусков алгоритма 35 запусков алгоритма 50 запусков алгоритма 
| т ке 
/ 
в = Знач. Знач. % г Знач. Знач. % —--- Знач. Знач. % 
ЭЛИТЫ улучш. ЭЛИТЫ улучш. ЭЛИТЫ улучш. 
28260 | 30936 | 33596 | 10.47 | 28198 | 30238 | 33623 27004 | 30246 | 33492 | 12.78 


33515 


Таблица 3 


особей 





Выбор элитной особи из пяти (граф — Ъга71158) 


Число 20 запусков алгоритма 35 запусков алгоритма 50 запусков алгоритма 
не ыы ыы 
/ 
а “_ Знач. | Знач. % —- Знач. | Знач. % —_ Знач. | Знач. % 
Й 
элиты | улучш. элиты | улучш. элиты | улучш. 
27635 | 30024 | 33264 | 12,76 | 27608 | 29885 | 33265 27667 | 29972 | 32948 | 11,72 


32745 


Анализ результатов вычислительного эксперимента позволяет построить диаграммы (рис. 1, 2, 3, 4, 5). 
14,66 


| 
а) 


6) в) 
Рис. 1. Процент улучшения элитной особи генетическим алгоритмом (выбор из трех). 
Количество особей и повторов: 750 (а), 250 (6), 500 (в) 





14,02 





В. Г. Кобак и др. Возможности использования элитных особей при решении задачи коммивояжера 


14,3 
| - В 
а) 6) в) 


Рис. 2. Процент улучшения элитной особи генетическим алгоритмом (выбор из пяти). 
Количество особей и повторов: 750 (а), 250 (6), 500 (в) 








20 


18,22 
17,84 


а) 6) в) 


Рис. 3. Процент отклонения среднего значения от оптимального (выбор элиты из трех). 





Количество особей и повторов: 750 (а), 250 (6), 500 (в) 
17,98 


ль 
15,85 


а) 6) в) 


Рис. 4. Процент отклонения среднего значения от оптимального (выбор элиты из пяти). 
Количество особей и повторов: 750 (а), 250 (6), 500 (в) 


6,4 








а) 6) в) 
Рис. 5. Процент отклонения лучшего значения от оптимального (выбор элиты из трех). 
Количество особей и повторов: 250 (а), 500 (6), 750 (в) 
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Анализируя результаты вычислительного эксперимента, можно отметить следующее: при увеличении количе- 
ства вершин в графе процент улучшения элитной особи увеличился, что говорит о повышении эффективности работы 
алгоритма с ростом числа вершин в графе. При выборе элитной особи из пяти с двукратным увеличением числа вер- 
шин в графе процент отклонения среднего значения от оптимума увеличился приблизительно на 6 %, что вполне до- 
пустимо. 

Для оценки эффективности работы алгоритма сравним его результаты с генетическим алгоритмом без элитной 
особи, но с той же стратегией отбора и таким же оператором мутации на графе Бга71 58. Результаты вычислительного 
эксперимента приведены в табл. 4. 

Таблица 4 
Результат работы генетического алгоритма без элитной особи 


Число повторов/ Среднее Лучшее % ОТКЛ. % ОТКЛ. 
особей из 35 запусков из 35 запусков среднего лучшего 


250 35223 29855 38,70 17,56 


500 33829 29827 33,21 17,45 
750 33716 28722 32,77 13,10 


Выводы. Ввод элитной муравьиной особи позволил уменьшить процент отклонения лучшего из найденных 





значений более чем в 2 раза, а процент отклонения среднего значения — приблизительно на 18 %. Увеличение числа п 
при выборе элитной особи позволило увеличить эффективность алгоритма примерно на 2 %. Следует также отметить, 
что с увеличением числа особей и итераций алгоритма улучшаются значения целевой функции. 
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